Resonant effects in the strongly driven phase-biased Cooper-pair box 
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We study the time-averaged upper level occupation probability in a strongly driven two level 
system, particularly its dependence on the driving amplitude xo and frequency u> and the energy 
level separation AE. In contrast to the case of weak driving (xo <C AE), when the positions of the 
resonances almost do not depend on xo, in the case of the strong diving (so ~ AE) their positions 
are strongly amplitude-dependent. We study these resonances in the concrete system - the strongly 
driven phase-biased Cooper-pair box, which is considered to be weakly coupled to the tank circuit. 

PACS numbers: 03.67.Lx, 03.75.Lm, 74.50.+r, 85.25.Am 



Several mesoscopic superconducting devices,which behave as quantum-mechanical two level systems (TLSs), were 
proposed and studied recently (see reviews ^, ^]). And although these devices are formally analogous to microscopic 
TLSs (such as electrons, atoms, photons, etc. [J]), they differ in that the coupling to controlling gates and an 
environment must be taken into account (this makes the numerical analysis of a mesoscopic TLS necessary). The 
study of the dynamic behaviour of the mesoscopic superconducting structures is interesting because they are suitable 
to observe the quantum-mechanical features by measuring macroscopic values and because of their relevance for 
engineering on the mesoscopic scale, e.g. for potentially realizable quantum computers based on superconducting 
Josephson o^ubits^The following non-stationary effects were studied in the superconducting effectively TLSs: Rabi 
oscillations 



qubits. The 



multiphoton excitations 



' euecis were 



12 . [13j , nonlinear excitations 



, Landau-Zener transition 

In this work we study the strongly driven superconductin g TL S. Namely, we study the phase-biased Cooper-pair 
box (PBCPB) (also called the Cooper-pair transistor) \lA Ilfil Il7l Il8l Il9| strongly driven via the gate electrode and 



probed by the classical resonant contour (tank circuit). The particular interest in this problem is because due to the 
interference between the Landau-Zener tunneling events, the system can be resonantly excited and the probability of 
the excitation quasi-oscillatory depends on the amplitude of the driving parameter 2(1 21 , 2^, . That is why we 
are interested in the dynamics of the strongly driven superconducting TLS - to clarify this problem and to relate it 
to the experimental results 2^| . 

The rest of the paper is organized as following. First we analyse the resonant excitations of a TLS, particularly, the 



difference between the weakly and strongly driven regimes. Then we study concrete situation of the strongly driven 
PBCPB, which is probed by the tank circuit. The paper ends with the conclusions. 
We consider a TLS described by the Hamiltonian 

H(t) = Aa x + (x f / + xq sinujt)a z . (1) 

Here a XiZ are the Pauli matrices. We are interested in the time-averaged upper level occupation probability, which 
is assumed to be related with the observable values. A driven TLS can be resonantly excited from the ground 



state to the upper state 



When the driving amplitude xo is small compared to the energy level separation 



AE = 2^J A 2 + x ^, the positions of the resonances in the time- averaged upper level occupation probability is 
determined by the multiphoton relation, AE — K ■ htu. Here to is the driving frequency and K is an integer. If the 
amplitude xq is increased, the position of the resonances is shifted (the Bloch-Siegert shift) Thus, at fixing u) 

and AE and with increasing the amplitude xo one should expect the (quasi-) periodical behaviour due to the shift of 
the multiphoton resonances. Below we analyse this issue in terms of the shift of the multiphoton resonances following 
Ref. 26] . Alternatively the quasi-periodical behaviour of the probability can be described in terms of the sequential 
Landau-Zener transitions with taking into account the quantum-mechanical interference between the transition events 



as in Ref. 



Consider first, for simplicity, the case of the zero offset, x ff = 0. In this case the position of the resonances in the 
dependence of the occupation probability on the system's parameters is defined by the following equation |26| : 

where E(k) = / dx ^ 1 Z^f^ - is the complete elliptic integral of the second kind and q — xq/A. The parameter q is 
J o v 

the key parameter of the problem; consider two limiting cases: of the weak driving, q <C 1, and of the very strong 
driving, q ^ 1 (the term "strong driving" we reserve for the case q ~ 1); from Eq. (J2J) it follows: 

AE = K ■ Tiuj, q < 1, (3) 
4rn 

-P- = ttK, <?» 1. (4) 
Too 

The first relation defines the multiphoton resonances, when the energy level separation, AE = 2 A , is a multiple 
of a photon energy htu. The resonances, defined by Eq. can be observed in the dependence of the occupation 
probability on lu or A, but not in the dependence on x . In the second case the resonances, defined by Eq. |0J, can 
be observed in the dependence on oj or xq, but not in the dependence on A; in this case equation Q also implies 



[2II and 3 



periodical (or quasi-periodical) dependence on the parameter <f> = Axo/fkv, which was studied in Refs. 
For the strong driving, q ~ 1, the resonances are expected in dependencies on each of the three parameters: to, 
xq, and A. Thus, we expect to find in the regime of the strong driving features typical for the two limiting cases: 
(i) quasi-periodic resonant dependence on xq and (ii) the resonances to appear in the dependence on A (with their 
positions being depende nt on X n 
Consider the PBCPB 
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16, 17 



23] excited through the gate electrode. The PBCPB is the small supercon- 
ducting island, which is connected via two Josephson junctions (characterized by energies -Eji,2 and phase differences 
61,2) to the ring with low inductance L (which is pierced by the magnetic flux <I> e ) and via the capacitance C g to the 
gate with voltage V g . The PBCPB is described by the Hamiltonian: 



,(0) _„(i) 



H=fd x -2E c (l-n g ' .,, 



sin ivt)a z 



(5) 



where the Coulomb energy of the island with the total capacitance Ctot is Ec = e 2 /2Ctot', the effective Josephson 

1 /2 

energy is e,j = [E^ + E 2 J2 + 2Ej\Ej2 cos <5) ; the total phase difference, 6 = 61 + 82, approximately equals to 
27r$ e /$o; and the dimensionless gate voltage is n g (t) — n g +n g sinwi = C g V g (t)/e. The Hamiltonian of the PBCPB 
(JSJ) coincides with the introduced above Hamiltonian with the substitutions: A = ej{S)/2, x ff = — 2Ec(l — n g ), 
and £0 = 2Ecn g ■ 

Now the parameter q is given by q = ^^-rig 1 ^. Thus both limiting cases - of weak and of very strong driving - 

n 

described above, can in principle be realized in the PBCPB [23J, where the domination of the Coulomb energy of a 
Cooper pair AEc over the coupling energy e,j is assumed, AEc/sj > 1. In Ref. ^| we have studied the case of weak 
driving and here we study the case of strong driving, q ~ 1, in detail. 

We will study the dependencies on and on 5 to demonstrate features (i) and (ii) . The occupation probabilities 
of the PBCPB are assumed to be probed by the tank circuit which is weakly coupled through the mutual inductance 
M to the PBCPB [27, 2|] . The average current (t) through the PBCPB is related to the phase shift between the 
voltage and current a, when the tank circuit with the capacity Ct and the inductunce Lt is driven at the resonant 



frequency lot = ^/VLtCt, as follows 



□ 



tan a ~ k 2 QL 



2e 



d(l 



h 85 



(G) 



where Q^ 1 = lotCtRt, k 2 = M 2 /(L ■ Lt)- To obtain the expectation value for the current in the qubit's ring, 
1^ = Sp (ijl^ ; we solve numerically the Bloch equations for the reduced density matrix p, as we did in Ref. nj . 
In Fig. ^ we plot the time-averaged upper level occupation probability P as a function of the amplitude n g x ^ at 



7r by making use of the solution of the Bloch equations. The case of ^ 1 (that is x a ff ^ 0) differs from 



the case of 



1 (x ff = 0) by the appearance of the additional peaks, which was discussed in Refs. 



2l| 



and 



J. We point out that similar dependence, which illustrates the feature (i), in the case x ff — can be calculated 
alternatively by making use of other approaches, namely with Eq. (13) from Ref. 2^] and with Eq. (17) from Rcf. 
[21 1. The numerical solution of the Bloch equations allows us to overcome the restrictions of the analytical works: in 
Ref. neither decoherence nor x a f / 7^ were taken into account while in Ref. | the assumption of very strong 
driving was done, which, for example, excludes the feature (ii). 




FIG. 1: Dependence of the time-averaged upper level occupation probability P (left) and phase shift a (right) on the amplitude 



n g at S — ir. 



Since at 6 — n the phase shift a is proportional to the time-averaged difference between the ground and excited 



state occupation probabilities 



3,1 



2P, Fig. ^presents also the dependence of a on n^. In Fig. |2]the dependence 



of the phase shift a on the total phase difference S is plotted for different amplitudes rig. Note that, as it was 
explained in Ref. the dependence of the phase shift a on S has hyperbolic-like character in the vicinity of the 

resonances. The parameters of the system taken for the Figs. Hand[5]are the following: Eji/Ec = 4.5, Ej2/Ec — 4, 
Tilo/Ec — 0.25, k 2 Q ■ 2e 2 LEc/h 2 = 0.01; the temperature was considered to be zero (i.e. much less than Ej\ — Eji)\ 
the decoherence rates we considered to be the functions of the energy level separation: r 12 oc AE(S) 1] (we have 
taken Tip/Ec ~ 0.01). Such set of parameters was taken for the comparison with the experimental results of Ref. 

y 

In conclusion, we clarified from analytical consideration the qualitative difference between the weak driving of a 
TLS and very strong driving. Then the strongly driven PBCPB was studied. The numerical results (Figs. n an d 
|2J demonstrated that (i) the dependence of the tank phase shift a on the amplitude n^P at 5 — ir has resonant 



quasi-periodic character and (ii) the resonances appear in the dependence on the phase difference S as the amplitude- 
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FIG. 2: The dependence of the phase shift a on the total phase difference 8 for different amplitudes . Upper curves are 
shifted vertically for clarity. 



dependent hyperbolic-like structures. We point out that the dependencies, characterized by the features (i) and (ii) 



similar to Figs. ^ and [3 were observed experimentally [2J|. And also similar to Fig. ^ quasi-periodic dependence 
of the upper level occupation probability on the driving (microwave) amplitude was observed in the superconducting 
TLS based on a large Josephson-junction in Fig. 6 of Ref. [5j- 

We thank the authors of Ref. ^| for communication of their experimental results prior to the publication and 
E. Il'ichev, W. Krech and V.I. Shnyrkov for fruitful discussions. The authors acknowledge the grant "Nanosystems, 
nanomaterials, and nanotechnology" of the National Academy of Sciences of Ukraine. 

Note added in proof. During the preparation of the manuscript we became aware that similar works on the 
strongly driven superconducting systems were done 2|1 3^|. The articles are given u p t o the experimental and 
theoretical study of the interference fringes in the strongly driven Cooper-pair transistor 
The conclusions of the present work agree with the conclusions of these articles. 



29j and the flux qubit 



anu 

Q. 
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